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Simple bimolecular reactions A\ + A2 f± A3 + A4 are analyzed within the framework of the 
Boltzmann equation in the initial stage of a chemical reaction with the system far from chemical 
equilibrium. The Chapman-Enskog methodology is applied to determine the coefficients of the ex- 
pansion of the distribution functions in terms of Sonine polynomials for peculiar molecular velocities. 
The results are applied to the reaction H2 + CI 5=^ HCl + H, and the influence of the non-Maxwellian 
distribution and of the activation-energy dependent reactive cross sections upon the forward and 
reverse reaction rate coefficients are discussed. 



I. INTRODUCTION 

Since the works of Prigogine and collaborators in the 1950's, the study of chemically reacting gases by means 
of the Boltzmann equation has constituted a topic of research [1, 2]. The first article in this series analyzed the 
chemical reaction A + A — > B + C under the assumption that the constituents have only translational energy and the 
reagents are more concentrated than the products [1]. The Boltzmann equation was solved by the Chapman-Enskog 
method, the distribution function was expanded to second order in the Sonine polynomials, and the reaction rate 
was determined for two kinds of reactive differential cross sections: one of them a step function, the other taking the 
activation energy into account . 

The second work analyzed the chemical reaction A^ + B ^ Ai + B and gave attention to the reaction heat [2]. The 
dependence of the reaction rate on the activation energy and on the reaction heat were determined. 

It has long been known that the reaction rate depends on the choice of the molecular interactions, and this motivated 
Present to propose the following expression for the reactive cross section [3, 4] 

0, E<e\ 

ttD 2 [1- (e*/E)}, E>e*. ' {l) 

where the formation of the activation complex is possible only when the distance between the molecular centers is 
equal to the parameter D, related to the diameters of the molecules. This cross section, known as the line- of- centers 
model, allows a chemical reaction whenever the relative translational energy of the molecules E is larger than the 
activation energy e*. 

In 1960 Ross and Mazur [5] analyzed the bimolecular chemical reaction A + B ^ C + D and extracted general 
expressions for the forward and reverse reaction rates and the entropy production from the Boltzmann equation. The 
distribution function was not determined explicitly from the Boltzmann equation so that no exact expressions for 
the reaction rates and entropy production were derived. The solution of the Boltzmann equation for the distribution 
function from the expansion of the distribution function in Sonine polynomials and computation of the forward and 
reverse reaction rate were presented later by Present [6] and by Shizgal and Karplus [7-10]. 

Many contributions to the analysis of chemically reactive systems follow these works. Among others we quote 
[11-32]. 

The cross sections determining the reactive collision term can be divided into two types, namely with and without 
activation energy [33, 34]. In general, cross sections with activation energy allow a reaction to occur whenever the 
relative translational energy exceeds the activation energy, as in Eq. (1). Even the slightest grazing collision leads to 
a chemical reaction. In a more realistic scenario, a reaction occurs only when the relative translational energy in the 
direction of the line joining the centers of the molecules is larger than the activation energy [11, 31, 34]. In this case 
the geometry of the collision plays a fundamental role in controlling the occurrence of a reaction. 

Another important issue in applications of the Chapman-Enskog to the reactive Boltzmann equation is the distortion 
of the Maxwellian distribution functions by the reaction heat [2, 14, 26, 31], which makes the reaction rate coefficients 
sensitive to the reaction heat. 
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The aim of this work is to analyze the influence of the activation-energy cross sections and of the non-Mawcllian 
distribution function on the reaction-rate coefficients. Although relying on the same method adopted in the above- 
mentioned works, we here solve the Boltzmann equation for the distribution function, which we expand in Soninc 
polynomials up to second order. The reactive differential cross section we consider allows reactions only when the 
relative translational energy in the direction of the line joining the centers of the molecules is greater than the 
activation energy, and we analyze slow reactions. In other words, we analyze the initial stage of a chemical reaction, 
in which the system is far from chemical equilibrium. In this stage the elastic collisions arc more frequent than reactive 
ones, and the affinity is much larger than the thermal energy of the mixture. We show that the reactive cross section 
under study markedly influences the reaction rates, while the effect of the non-Maxwellian distribution is not too 
pronounced. 

The work is structured as follows: Section 2 introduces the system of Boltzmann equations for a simple bimolecular 
chemical reaction A\ + A 2 ^ A3+A4. The expressions for the reaction rate coefficients that follow from the Boltzmann 
equation and the specification of the elastic and reactive cross sections are the subject of Section 3. In Section 4 
the Arrhenius equation is obtained from a Maxwellian distribution function and it is shown that the reaction-rate 
coefficient when the relative translational energy in the direction of the line which joins the centers of the molecules 
must be larger than the activation energy is different from the coefficient when only the relative translational energy 
must be greater than the activation energy. The analysis of the slow chemical reactions is the subject of Section 5, 
where the distribution functions are expanded in Sonine polynomials of the peculiar molecular velocity and the first 
coefficients of the expansions are obtained with the Chapman-Enskog methodology. In Section 6 the results of the 
previous section are applied to the chemical reaction i? 2 + CI f± H CI + H, and the coefficients of the non-Maxwellian 
distribution function, and the forward and reverse reaction coefficients are plotted as functions of the temperature. 
The results show the effect of the non-Maxwellian distribution functions and of the reactive cross sections on the 
reaction-rate coefficients and on the entropy production rate. 



II. BOLTZMANN EQUATIONS 

We consider a simple reversible bimolecular gas reaction characterized by the chemical law A\ + A 2 ^ A3 + A4, 
which takes elastic and reactive binary encounters between the molecules into account. 

The elastic collisions between the two constituent molecules, of masses m a and mp, have asymptotic pre-collisional 
velocities (c Q , c^), asymptotic post-collisional velocities (c^, c^), and asymptotic relative velocities gp a = cp — c a and 
g,'p a = c'p — c' a , respectively, so that the conservation laws of linear momentum and energy are given by the expressions 

m a c a + mpep = m a c' a + mpc'p, (2) 

1 2,1 2 1 '2,1 12 (o\ 

-m a c a + -mpcp = -m a c a + -mpep , (3) 

respectively, where a, (3 = 1,...,4. The energy conservation law can also be written in terms of the asymptotic 
relative velocities as gp a = g'p a . 

For a reactive collision the conservation laws of mass, linear momentum and total energy read 

mi + rn 2 = m 3 + m 4 , (4) 
miCi + m 2 c 2 = ?n 3 C3 + m 4 c 4 , (5) 

ei + \fnic\ + e 2 + ^m 2 cj = e 3 + \^m z c\ + e 4 + ]^m A c 2 A , (6) 

respectively. Here e a denotes the molecular binding energy (a = 1,...,4), while (04,02) are the velocities of the 
reactants, and (C3, C4), the velocities of the products of the forward reaction. The conservation law for the total energy 
(6) can be written in terms of the relative velocities gp a = \gp a \ = \cp — c a \ and of the heat of reaction Q — defined 
as the difference of the chemical binding energies of the products and the reactants Q = e 3 + e 4 — e x — e 2 — as 

2-TO12321 = 7^34343 + <2> ( 7 ) 

where m a p = m a mp/ (m a + mp) denotes the reduced mass. 

We characterize the state of a reacting gaseous mixture in the phase space spanned by the positions x and velocities 
c a of the molecules by the set of distribution functions f a = /(x, c Q ,i) with a = 1, ... ,4. The distribution function 
f a is defined so that the number of a molecules in the volume element dx.dc a around the position x and the velocity 
c Q at time t is given by f a dxdc a . 
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The phase-space evolution of the distribution function f a for constituent a is governed by the Boltzmann equation, 
which in the absence of external forces reads 



+ c fjh = E / ti'J'p - UU) 9p^ a dO dc /3 + Q« a = 1, . . . , 4, (8) 
0X1 0=x J 



dfa a df, 
dt 



with the shorthand f' a = f a (x,c' a ,t). 

The left-hand side of Eq. (8) refers to the space-time evolution of the distribution function, while its right-hand side 
takes the molecular collisions into account. The latter has two terms. The first one describes the elastic interactions 
among the constituents. In the integrand, the factor a a p is the differential elastic cross section, and df2 = sm\dxde 
the solid-angle element, with x denoting the scattering angle, and e, the azimuthal angle that characterizes the 
collision. 

The second term on the right-hand side of Eq. (8) is related to the reactive collisions. The expression for the 
constituent labeled by the index 1 is obtained as follows. The number of reactive collisions for the forward reaction 
Ai + A2 — + A 4 per unit of volume and time is given by (/1 f 2 (721 c*2 dQ,* dc% dc 2 ), where a\ 2 represents the 
reactive differential cross section and dfi* is the solid-angle element relating the orientation of the post-collisional 
relative velocity g43 with respect to the pre-collisional one g2i. Likewise, the number of collisions for the reverse 
reaction A\ + A 2 t— A3 + A 4 reads (/3 f 4 g 43 CT34 dtt* dc 3 dc 4 ). For a fixed value of the heat of reaction it follows 
from Eq. (7) that m\ 2 g 2 \ d#2i = 77734(743 d<?43, and we also may obtain the equality 771341721 dc3 dc 4 = m\ 2 g 43 dci dc2. 
Furthermore, for small heats of reaction the principle of microscopic reversibility yields (77712(721 ) 2 Ci 2 = (77734343) 2 C34. 
Hence, all relationships accounted for, we find that f 3 f 4 343 1734 dft* dc$ dc 4 — f 4 (77712/77734) 3 321 cr\ 2 ^* dci dc2. 
From the difference of the expressions for the number of collisions in the backward and forward reactions divided by 
dci and integrated over all values of dc2 and df2*, we obtain the following expression for the reactive collision term 
for the constituent labeled 1: 



\m 34 J 



a* 12 g 2 idn*dc 2{1) . (9) 



The reactive collision term for the constituent labeled by the index 2 being similar, the terms within parentheses in 
Eq. (9) determine Q 2 . In the same notation, the reactive collision terms for the constituents labeled by the indexes 
3 and 4 read 



qR 

^3(4) 



/1/2 ( — — ] - /3/4 



CT34343dfi*dc 4 ( 3 ) . (10) 



III. REACTION RATES AND DIFFERENTIAL CROSS SECTIONS 

The description of the evolution equation for the particle number density n a — J fa dc a of constituent a is obtained 
by integrating the Boltzmann equation (8) over all values of c Q , which yields 

<977 a ^ dn a vf _ ^ 
dt dxi a 

Here, the bulk velocity vf and the particle number density production r a of the constituent a read 

V" = / Cf f a dc a , T a = -V a (773774k,. - 77 x 772k/) , (12) 

77ck J 

where kf and k r denote the forward and reverse reaction rate coefficients, respectively. They arc defined by 

k f = — !— [ fif20-* 2 92idQ*dc 1 dc 2 , k r = ( hh ( — J ^ 2 32idO*dcidc 2 . (13) 

77i 77 2 J n 3 n 4 J \m 34 J 

In Eq. (12)2 we have introduced the stoichiometric coefficients v a of the constituent a, which for the chemical reaction 
Ai + A 2 ^ A 3 + A 4 are given by v\ = v 2 = — v 3 = — v 4 = — 1. 

In order to determine the distribution functions f a from the system of Boltzmann equations (8), we have to specify 
the elastic and the reactive differential cross sections. We assume that the elastic differential cross sections correspond 
to a hard-sphere potential, i. e., that 



iip, d Q(3 = - (d Q + d^) , (14) 
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where d Q and represent the diameters of the colliding spheres. 

For the reactive differential cross section we shall use the modified line-of-centers model [If, 31, 34] 



0. 



'34 



is H 2 



2e f 



1 



mi2(k 2 i-g2l) 



2e r 



TO34(k43-g43) 



"ii2 (k 2 i • g2i) /2 < e/, 

, "112 (k 2 l • g2l) 2 /2 > £/, 

™34 (k 4 3 • g43) 2 /2 < e r , 
, »™34 (k43 ' g43) 2 /2 > £ r . 



(15a) 
(15b) 



Here e/ and e r denote the forward and reverse activation energies, respectively, and Sf and s r the corresponding steric 
factors, while k 2 i and k43 are unit collision vectors in the directions of the centers of the colliding molecules pointing 
from the centers of 2 and 4 to the centers of 1 and 3, respectively. 

The differential cross sections in Eqs. (15a) and (15b) take into account the activation energies and the geometry 
of the reactive collisions, since they arc functions of the relative translational energies in the direction of the lines 
joining the centers of the molecules ran (k 2 i • g 2 i) 2 and 77134 (k43 • g43) 2 . The original line-of-centers model considers 
the relative translational energies as m,\292\ an< ^ to 34<?43j so that even the slightest grazing collision leads to a chemical 
reaction. 



IV. ARRHENIUS EQUATION 



For mixtures the chemical potential of the constituent a, which takes into account the binding energy of the 
molecules — but not the internal states of the molecules associated with the rotational, vibrational, electronic and 
nuclei states — is given by the equality 



Ma 



kT 



|^ + In n a - ^ In T + C a 



(16) 



where k denotes Boltzmann's constant, C a is a constant, and all constituents have been assumed to be at the tem- 
perature T of the mixture, defined by 



T = — — [ m a \c a - v a \ 2 f a dc a . 
3kn a J 

Chemical equilibrium is characterized by the condition 

4 

^ rn a u a ^ = 0, 



(17) 



(18) 



where the index " eq" denotes the equilibrium- value of the chemical potential. Out-of-equilibrium chemical reactions 
are described by the affinity, defined as 



A = — ^] rn a v a ^, a , with _4 cq = 0. 



(19) 



a=l 



From the above equations we can obtain the law of mass action, 



In 



"-1 "-2 

cq cq 1 
n 3 n A \ m 



3 

/ m 3 m 4 \ 2 
im 2 J 



and the following expression for the affinity 



In 



cq cq 
ca cq" 



A\ 



(20) 



(21) 



where we have introduced the dimensionless reaction heat Q* = 
energy kT. 



Q/kT and affinity A* = A/kT, in units of the thermal 
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The state of equilibrium of a non-reacting gas mixture is characterized by the Maxwcllian distribution functions 



•' O: 



(0) 



/ m a 



V 



a \2vkT) 



a 



(22) 



where the particle number densities are uncorrelated. 

In fact, these distribution functions make the elastic collision terms of the Boltzmann equation (8) equal to zero. 
In general, however, the reactive collision terms (9) and (10) do not vanish, because the particle number densities in 
the distribution functions (22) are not the equilibrium densities, which characterize the chemical equilibrium and are 
related by the law of mass action, Eq. (20). In Eq. (22) we have introduced the peculiar velocity £f = cf — Ui where 

i'i = Yltt=i Qa v i ' I ELi Qa is the bulk velocity of the mixture. 

From the Maxwellian distribution functions we can calculate the forward reaction rate coefficient by substituting 
Eq. (22) in the definition of the forward reaction rate coefficient, Eq. (13). Substitution of the differential cross 
section (15a) and integration of the resulting equation leads to the result 



,(°) 



(mim 2 y- 



exp 



(mi 



+ m 2 )G\ 2 _ mi 2 gli 
2kT 2kT 



Sfdi 2 



2c, 



mi2 (.921 cos6>)^ 



,g2i sin # cos # d0 de dg2i dGi 2 , (23) 



since the solid-angle element is dO* = 4 sin 9 cos 9 d9 de where 9 and e are the angles characterizing the reactive 
collision process, while k 2 i • g 2 i = <?2icos0. We have introduced the relative velocity g 2 i and the center of mass 
velocity G12 defined by the relation 



a/3 



Si Si 1 



G 



a/3 



(24) 



The range of the integrals in the variables G\ 2 and e on the right-hand side of Eq. (23) are < G\ 2 < 00 and < 
e < 2tt, respectively. The range of the integral over the angle 9 follows from Eq. (15a), according to which the normal 
relative velocity has a lower bound so that the interval of integration of 9 is given by < 9 < arccos ^/2e//mi 2 g| 1 . 
The range of the relative velocity is y / 2e//mi2 < 512 < 00. 

From the integration of (23) we obtain the following expression for the forward reaction rate coefficient 
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l 8nkT sjdl 2 e.-tf 



(25) 



a modified Arrhenius equation. Here, is the forward activation energy in units of the thermal energy fcT, and 
Ei(eJ) represents the exponential integral Ei(e^) = f°? e V y dy . 

The underlined term in Eq. (25) is the usual Arrhenius equation, which follows from substituting mi 2 g 21 for 
WI12 (k2i • g2i) in the differential cross section. Inspection of Eq. (25) shows that the reaction-rate coefficient is 
smaller for the latter case, in which a reaction occurs only when the relative translational energy in the direction of 
the line which joins the molecular centers exceeds the activation energy. In the former case even a grazing collision 
with relative translational energy larger than the activation energy would lead to a chemical reaction. 



V. SLOW CHEMICAL REACTIONS 



When dealing with chemical reactions within the framework of the Boltzmann equation, we can analyze two regimes: 

Fast reactions: the last stage of a chemical reaction. The system is nearby a chemical equilibrium state. The 
frequencies of the elastic and reacting collisions are of the same order of magnitude and the affinity is small in 
comparison with the thermal energy of the mixture, i.e., A <C kT; 

Slow reactions: the initial stage of a chemical reaction. The system is far from chemical equilibrium. The clastic 
collisions are much more frequent than the reactive ones and the affinity is much larger than the thermal energy 
of the mixture, i.e., A ^> kT. 
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In this work wc shall analyze slow reactions. We write the Boltzmann equations in the form 



Vf a + A£ 



dfa 

dxi 



OXi 



hh 
hh 



I mi2 
\ m 12 



- hh 



hh 



l 4 f 

a* 2 g 21 (Kl*dc 7 = / {f' a f'p - fafp) gp a <rp a d£ldcp, (26a) 

0=1 J 

1 4 f 

o-| 4 g 43 dr2*dc 7 = -^] / (/4/g - / a /a) gpaVpadttdcp, (26b) 



where 2? = <9/di + Vid/dxi is the material time derivative, and A a parameter of the order of the Knudscn number. 
In Eq. (26a) a = 1, 2 corresponds to 7 = 2, 1, while in Eq. (26b) a = 3, 4, to 7 = 4, 3. Since we assumed the elastic 
collisions to be more frequent than the reactive ones, Eqs. (26a) and (26b) indicate that the reactive collision terms 
and the material time derivatives are of the same order of magnitude, while the gradients of the particle number 
densities, velocity and temperature of the mixture are of the next order. 

To obtain the distribution functions from Eqs. (26a) and (26b) we apply the Chapman-Enskog methodology. To 
that end, we expand the material time derivative in a power series of a parameter A, i. e., we write 



V = V a + \T>i + A 2 £> 2 

along with the distribution functions 

f a = fL 0) + xfL 1) + x 2 f a 2) + -.., 



a = 1, 



.4. 



(27) 



(28) 



The parameter A in the distribution function will be set equal to the unit later. 

We insert the expansions (27) and (28) into the Boltzmann equations (26a) and (26b). Comparison of equal powers 
of A leads to the following system of integral equations: 



o = E / (f^fP'-fL 0) fP)s^ Pa dndc,, 
0=1 J 

/r/ ( ) (^) 3 -/f)/H a* 12g21 dn*dc 7 =i: f + 

\ ''"34 / B—l 



(29) 



f^f^-f^f^Qp^padSldcp, a =1,2 7 = 2,1, 



(30) 



V f^ - 



a (0) / 2 (0) (2*) 3 - /i 0) /i 0) ^dwdc, = £ / [/M'/r + /r /r 



-fPfr-fZ'fngpaVpadtodcp, a = 3,4 7 = 4,3. 

(31) 

Here we have only considered the integral equations to first order in A, i. e., associated with the zeroth order term 
/a°^ and the first-order correction . 

The solution of the integral equation (29) is the Maxwellian distribution function (22). To solve the integral 

equations (30) and (31) we assume that the reaction heat affects the Maxwellian distribution functions, so that fa 
can be expressed in terms of Soninc polynomials in the peculiar velocity i. e., 



J a — J a / < dna °1 



(n) ( m a £,l 
2kT 



where Si (x) = 



T(n + 3/2) 



fe=0 



k\(n-ky.T(k + 3/2) 



(-*)* 



(32) 



The coefficients a nQ are constants to be determined from the integral equations. Up to second order terms in the 
expansion, the distribution function can be written as 



Ja — J a 



(0) 



1 + ai 



3 _ m a £, a 
2 2kT 



a 2a 



15 5m a ^ m a £ a 



AkT 



8(*T) 5 



a = 1, . 



(33) 



Insertion of the distribution function (33) into the definition (17) of the temperature of constituent a, shows that 
all coefficients ai Q are zero. To determine the remaining coefficients a2 Q we multiply Eqs. (30) and (31) by the Soninc 
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polynomial [15/8 - hm a £ 2 JAkT + m 2 ^/8(fcT) 2 ] and integrate over all dc a , which shows that 



1-e" 



s f 2 di 2 e- e > l-4e). + 3e* f Ei(e* f )e e > 



m{ 2 nin2 



8ml 



n a np^/rn^ ^ ( 10m 2 a + 8m a mp + 13to^ ^ 
^(rria+mp) 2 a/3 | m a +711/3 



a = 1,2; 



;i - e-- 4 * ) sfd^e" 6 } {l - 4 e) + 3 e) Ex (e)) e 6 / - 4Q* [l + Q* - e} Ei (e)) e e / (3 + Q*)l } mj 4 rnn 2 



8m 2 m £ 2 
10m 2 + 8m a mp 



13m 2 



^ (m a + mp) 



2 



Too, + m/3 



a 2 Q - 15m a/ 3a 2/ 3 



(34) 



3,4. (35) 



From the algebraic system of Eqs. (34) and (35) we can determine the four coefficients a 2 i, . . . , a 2 4 of the distribution 
functions (33). The expressions for these coefficients are too long to be presented here. 

To determine the forward and reverse reaction rate coefficients, we insert each distribution function (32) in the 
corresponding definition, in Eq. (13). Integration of the resulting equalities then yields the expressions 



</ = k? 1 



i-4 t ; + 3<;E,( t ;)e'; m«j-fnnj 



l-e* f Ei(e})e e / 



mi + m 2 ) 2 



(36) 



(nwa\* 0%( o) J 1 1 ~ 4 ( C ) + Q* + Q* 2 j + g Me}) e'/ (3 + 12Q* + 4Q* 2 ) ^ + ^2 

£ } El (e}) 



Vm 3 m 4 y J I 1 - e *Ei(e*,)e e / 8(m 3 + m 4 ) 2 



The underlined terms are corrections to the forward and reverse reaction rate coefficients, for non-Maxwellian 
distribution functions. The corrections depend on the coefficients a 2 i, . . . , a 2 4 of the distribution functions (33), on 
the reaction heat, and on the activation energy of the forward reaction. 



VI. REACTION RATES FOR H 2 + Cl^ HCl + H 



We now apply the results of Sees. IV and V to evaluate the rate coefficients for the reaction H 2 + CI HCl + H. 
To that end we need to know certain characteristic parameters for the constituents of the mixture, such as masses, 
diameters, forward and reverse activation energies, reaction heat and steric factors. Table 1 lists the molecular weights 
M a and the coefficients of shear viscosity fx a at temperature T = 293 K for the single constituents H, H 2 , CI and 
HCl [35]. From the expression of the coefficient of shear viscosity for hard-sphere molecules (sec e.g. [34]) 

^ = l6djV^' (38) 

we can calculate the diameters of the single constituents H 2 and HCl. For the constituent CI we take the diameter 
to be twice the atomic radius, whereas for the constituent H we take the diameter to be of order of twice the Bohr 
radius a Q = 0.529 x 10~ 10 m. 



Gas 


H 


H 2 


CI 


HCl 


M a 


1.008 


2.016 


35.453 


36.461 


Ho, (xl0~ 5 Pa s) 




0.841 




1.332 


d a (xltT 10 m) 


1.06 


2.78 


1.99 


4.55 



TABLE I. Molecular weights, viscosities and molecular diameters. 
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Tabic 2 lists the coefficient A of the Arrhenius equation k^ ' = Ae~ e f, the forward, and the reverse activation 
energies [36]. The reference temperature for this reaction is 300 K and the reaction heat was obtained from the 
relation Q = e/ — e r . The steric factors in this table were calculated from the identification A = yJ%-xkT jm\2 sf d^ 2 , 
from Eq. (25). 



Reaction 


A(m 3 /mol s) 


e/(kJ/mol) 


e r (k J /mol) 


Q(kJ/mol) 


s / 


H 2 + CI ->• HCl + H 


7.94X10 7 


23.03 


18.84 


4.19 


0.648 


HCl + H ->■ H 2 +Cl 


4.68 xlO 7 


18.84 


23.03 


-4.19 


0.497 



TABLE II. Arrhenius coefficients, forward and reverse activation energies, reaction heat and steric factors. 



To analyze a slow reaction, in which the concentrations of the reagents is larger than those of the products, we start 

out with a preliminary evaluation of the factor appearing in the reverse reaction rate coefficient (37). 

If we adopt the tabulated values and consider the reaction in the direction H 2 + CI HCl + H, a very large ratio 
results between the reverse and the forward reaction coefficients. We therefore discard this alternative and consider 
only the reaction in the direction HCl + H — * H 2 + CI, for which the ratio is small. 

To determine the coefficients a2 Q from the algebraic system of Eqs. (34) and (35) we need to know the temperature 
and the concentrations of the constituents x a = n a /n where n = X^=i Tip and ^* =1 x a = 1. Here we analyze the 
case x\ = X2 and X3 = x<±. Figure 1 shows the coefficients a2 a , plotted as functions of the temperature in the range 
400 K < T < 1200 K for two concentrations: X\ = 0.30 and X\ = 0.35. The figure shows that the deviations from 
the Maxwellian distribution functions — which are given by the coefficients a2a — increase with the temperature. In 
comparison with the departure for the constituent H 2 , the departures for the constituents HCl, H and CI are not 
very large. 

Figure 2 shows the dependence of the reaction rate coefficients on the temperature in the range 400 K < T < 1200 K, 
for two concentrations: X\ = 0.30 and X\ = 0.35. The left (right) panel represents the forward reaction (reverse 
reaction) rate coefficient. Consider first the influence of the collision geometry on the reaction rate coefficients. The 
dashed-dotted lines in the inset left panel and in the right panel represent the Arrhenius law, given by the underlined 
term in Eq. (25), while the solid lines correspond to the modified Arrhenius equation (25), when the geometry of the 
collision is taken into account, i.e., when a reaction occurs only if the relative translational energy in the direction of 
the line joining the molecular centers exceeds the activation energy. The reaction rates in the former case are larger 
than those in the latter, an indication that more reactions occur in the former case, due to the fact than even grazing 
collisions with translational energy larger than the activation energy lead to chemical reactions. In the same figures we 
also plot the influence of the departures from the Maxwellian distribution functions on the reaction rate coefficients. 
The left and the right panels show that the non-Maxwellian distribution functions reduce the reaction rate coefficients 
and that the forward reaction rate coefficient is less reduced than the reverse reaction one. The reduction becomes 
more pronounced as the reagent concentration grows. 

Once the forward and reverse reaction rates are known, we can compute the entropy production rate. According 
to a phenomenological theory (see e.g. [37, 38]) the entropy-density rate <r is given by 

■; = — — , where — = k f nin 2 - k r n 3 ri4, (39) 

the so-called reaction velocity. 

Equation (39) can also be deduced from the Boltzmann equation; for more details see Ref. 5 and 34. In view of 
the law of mass action kfn^n^ = k r ng q n^ q , which shows that, in equilibrium, the forward and the reverse reactions 
occur with the same probability, the affinity (21) becomes 

A = kT\n( W J^\, (40) 
V k r n 3 n 4 I 



so that the entropy production rate (39) can be written as 

J_ = fcl nf^ 



h = fcln ( rSS ) {k f XlX2 - KxzXa) ■ (41) 
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FIG. 1. Coefficients a2 Q as functions of the temperature for two concentrations: 0.30 and 0.35. 




As one would expect, the entropy production rate is positive semi-definite — thanks to the inequality (x — 1) In a; > 
which is valid for all x > 0, the equality being valid for x = 1. 

Figure 3 shows the entropy production rate (41) as a function of the temperature for the Arrhenius and modified 
Arrhcnius cases, to show that in the former case the rate is larger than in the latter. This is expected, given the 
entropy production rate dependence on the reaction rates. For larger concentrations x\ the entropy production rate 
increases in both cases. 



VII. CONCLUSIONS 



We have used the Boltzmann equation to analyze a bimolecular chemical reaction in its initial stage. As an 
illustration, we computed the reaction rates and the entropy production rate for the typical bimolecular reaction 
H-2 + CI f± HCl + H . To solve the Boltzmann equation by the Chapman-Enskog method we have expanded the 
distribution function in Sonine polynomials up to the second order. Our reactive differential cross section allows a 
reaction to occur only when the relative translational energy in the direction of the line of centers which joins the 
centers of the molecules is greater than the activation energy. With this cross section the calculated reaction rates 
are substantially smaller than those obtained from the line-of-centers model, which requires the relative translational 
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FIG. 3. Entropy production rate in J dm 3 



/(K mol s) as function of the temperature in K. 



energy to be larger than the activation energy. We showed that the restriction to large relative translational energies 
along the line of centers affects the calculated rates more than the departure from the Maxwcllian distribution function. 
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